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Abstract. Let {Ge)!iyo be a family of 'e-thin' Riemannian manifolds modeled 
on a finite metric graph G, for example, the £-neighborhood of an embedding 
^ 1 of G in some Euclidean space with straight edges. We study the asymptotic 

Q ' behavior of the spectrum of the Laplace-Beltrami operator on Ge as e ^ 0, for 

' various boundary conditions. We obtain complete asymptotic expansions for 

, the fcth eigenvalue and the eigenfunctions, uniformly for k < C£~^ , in terms of 

■ scattering data on a non-compact limit space. We then use this to determine 

the quantum graph which is to be regarded as the limit object, in a spectral 
sense, of the family (Ge). 

Our method is a direct construction of approximate eigenfunctions from 
. the scattering and graph data, and use of a priori estimates to show that all 

' eigenfunctions are obtained in this way. 
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1. Introduction 

Consider a graph G with a finite number of vertices and edges embedded in R" 
with straight edges. For e > let be the set of points of distance at most e from 
G. For small e, Gg 'looks almost like' G; in other words, the one-dimensional object 
G should be considered as a good model for G^, and so one expects many physical 
and analytical properties of G^ to be understandable, in an approximate sense, 
by corresponding properties of G. The property we analyze in this light here is 
the spectrum of the Laplacian, under various boundary conditions. For motivation 
from physics, see for example |13j . 

This problem has received some attention in the last decade. The case of Neu- 
mann boundary conditions was analyzed at various levels of generality in [4], [6], 
[Is] . [21], [5], where it was shown that for each fc e N the /cth eigenvalue of the 
Neumann Laplacian on G^ converges, as e — > 0, to the /cth eigenvalue of the second 
derivative operator on the union of the edges of G, where at the vertices so-called 
Kirchhoff boundary conditions are imposed. The question of what the correspond- 
ing limiting behavior is for other, for example Dirichlet, boundary conditions was 
characterized as 'very difficult' in 5 and remained open until some partial progress 
was made recently, see Section [1^41 

In this paper we solve this problem for a general mixed boundary problem, 
where we impose Dirichlet boundary conditions on one part of the boundary and 
Neumann conditions on the rest. Other boundary conditions, for example Robin, 
are also possible. Instead of the setting of an embedded graph described above, 
we consider the more general, and mathematically more natural, situation of a 
shrinking family of Riemannian manifolds modeled on the graph: We now consider 
G as an abstract metric graph, that is, for each edge e a positive number If, (to 
be thought of as half the edge length) is given. In addition, we have geometric 
data: For each edge e an (rt — l)-dimensional Riemannian manifold YJ,, and for 
each edge v an n-dimensional Riemannian manifold X^] all these manifolds are 
compact and may have a piecewise smooth boundary; also, for each edge e incident 
to a vertex v (denoted e ~ w), we are given an identification (gluing map) of Ye 
with a subset of the boundary of X^, without overlaps. Then Gg is defined by 
gluing cylinders of length 2lf. and cross section eYg (the factor denotes a rescaling 
of the metric) to eXy for all pairs e ^ v. Also given are subsets D and N of the 
boundaries of each X^ and Yg, which yield corresponding subsets of Gg (or, for 
Robin conditions, similar data). We investigate the Laplace-Beltrami operator on 
Ge , where Dirichlet boundary conditions are imposed on D and Neumann conditions 
on N. We suppose sufficient regularity of the boundary and the D/N decomposition 
so that this problem has discrete spectrum. 

Besides the e-neighborhoods mentioned before this manifold setting includes the 
case of the boundary of the e-neighborhood of G (in this case, G^ itself has no 
boundary). Our results imply the following theorem. 



Theorem 1. Let ^ii{e) < M2(£) !i • • ■ be the eigenvalues of the Laplacian on Ge, 
with the given boundary conditions, repeated according to their multiplicity. 

For each edge e, let Ve be the smallest eigenvalue of —A on Yg, with the given 
boundary conditions, and let v = miUg v^, . 
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There are numbers D Cz No, < ti < T2 < ■ ■ ■ < td < v and < &i < ^2 < 
• 00 so that for e — > 

(1) A^fe(e) =e-Vfc + 0(e-^/^), k^l,...,D 

(2) ^k{e)^e-^iy + bk-D + 0{e), k>D 

In fact, we obtain much more precise information: complete asymptotics up to 
errors of the form 0{e^'^^'^) also in the case k > D, uniform estimates not only for 
fixed k but also for k on the order of as well as detailed information on the 
eigenf unctions, see the theorems below. 

For Theorem [1] to be meaningful it is essential to identify the numbers Tk and bk 
in terms of the given data. We do this in two steps: First, we give general formulas 
identifying these numbers in terms of scattering data of an associated non-compact 
n-dimensional limiting problem. Second, we analyze how this scattering data can 
be obtained from the given data. It turns out that in this second step the situation 
of general boundary conditions carries some essential new features compared to 
Neumann boundary conditions: while in this case one has D — and the bk are 
determined by the graph and the edge lengths alone (that is, independent of the 
manifolds Xy, Yg), in general D > and the and bk depend also on transcendental 
analytic information involving the manifolds Xy and Ye'. The are L^-eigenvalues 
in the afore-mentioned limit problem, and the bk are eigenvalues of a quantum 
graph whose vertex boundary conditions involve the scattering matrix of that limit 
problem, see Theorems [2] and [3l 

1.1. Main results. We now state our results more precisely. Because of the diver- 
gence in ([T]), ([2]) it is more convenient, and for a proper understanding it turns out 
to be essential, to rescale the problem: We multiply all lengths by e~^. We denote 

N:=e-\ X^:=e-'G, 

This rescales the eigenvalues by the factor e^. In Xq , the vertex manifolds and 
the cross sections of the edges are independent of N while the lengths of the edges 
are 2Nle, and we are interested in the limit N 00. Central to the analysis is the 
limit object, to be thought of as limAr^oo j 

(3) X°° ■.= [JX^ (disjoint union) 

V 

where the 'star' X!^ of a vertex u of G is obtained by attaching a half infinite 
cylinder [0, 00) x to the vertex manifold Xy, for each edge e incident to v, see 
Figure [H Denote by Y the cross-section for X°°, i.e. the disjoint union of the Ye, 
each one appearing twice (once for each endpoint). The graph structure is encoded 
by a map a : Y ~* Y which toggles the two copies of Yg for each edge e. See Section 
[2] for precise definitions. 

Since X°° is a non-compact space, the spectrum of the Laplacian — Ax<» (with 
corresponding boundary conditions) is no longer discrete. Its spectral theory is 
obtained via scattering theory, that is, one considers it as compact perturbation of 
[0, 00) X Y, the disjoint union of the edge cylinders. Let v < vi he the smallest 
and second smallest eigenvalue of the Laplacian of Y and let V be the eigenspace 
corresponding to v. The absolutely continuous spectrum of — Ax°° is [1^,00). To 
each \ = V + ^ [v,vi) IS, associated the scattering matrix S{a), a linear map 
V ^ V depending holomorphically on a, and to each Lp £V a. scattering solution 
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(generalized eigenfunction) E^^^. These are the main objects of scattering theory. 
In addition, — Ajcoo may have discrete spectrum (L^-eigenvalues). 

We also have linear maps L,a : V ~^ V where L encodes the edge lengths (it 
equals /gld on the subspace of V corresponding to the edge e) and a is induced 
by the graph structure map above. The maps a and 5(0) are involutions on V, so 
they have only the eigenvalues ±1. 

We need the following data derived from the spectral data of — Ajfoo: 

• The L^-eigenvalues t[ < T'2 < . . . < t'^ < vi —Ax-^, with eigenspaces 
fi,...,£rfCL2(X°°); 

• the space Vq Ker(Id - cr) n Ker(Id - 3(0)) C V. 

• the spaces Vz '■— Ker(Id — e^^^'"aS{0)) C V for z > 0. This is non-zero for 
a discrete set 

(4) 0<zi<Z2<-.- the zeroes of /(z) = dct(Id- e'^2L^5.(-0)) 

• for each i € N, pairwise different holomorphic functions Z^{a), p = 1, . . . , 
with Ti g N, defined on C \ {a S K. : |a| > vi} and with 

(5) (0) - Zr 

for each p, and holomorphically varying subspaces C V with 0^'^^ '^ii^) 

The and arise from a bifurcation analysis of Ker(Id — e^^'"^^ a S {a)) , see 
Theorem [121 

Theorem 2 (Main Theorem). Denote by —Axn the Laplace- Beltrami operator 
on Xq , with the given boundary conditions. For any Amax < t^i there are c,No>Q 
such that for N > Nq the spectrum of — Aj^-n in [0, Amax] consists of (counting 
multiplicities ) 

a) dimfp many eigenvalues of the form Tp + 0(e~'^^), for p — 1, . . . ,d, 

b) dimT'o many eigenvalues of the form v + 0{e~''^), 

c) diniT'f many eigenvalues of the form 

(6) + ™i;,/3^1zf(l) + 0(e-^) 
for each i e N and p — 1, . . . , for which v + (3"^ < Amax- 
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The constants c, Nq and the constants in the big O only depend on the spectral data 
of Ax°° ■ In particular, the asymptotics are uniform for eigenvalues < Amax- 

In terms of e = A^^^ ([6|) gives an eigenvalue fi = u + e^zf + 0{e^) by ([5|), so 
this implies Theorem [T] take D — X]p=i*i™'^p with d' = max{p '■ t'^ < v}\ the 
Tfe are the repeated dimfp times; and the bk are repeated dimT'o times, and 
the numbers zf repeated A\mVi many times. The Tp with p > d' are invisible in 
Theorem [1] since there k is fixed as e ^ 0. 

Theorem [5] also describes eigenvalue fik with /c ^ oo as £ — > 0. We determine 
the permitted range of k. From /x^ = j/ + e'^zf + 0{e^) < Amax one obtains ezi < 
V Amax — 1^ + 0{e), and the Weyl asymptotics for the Zi, see Proposition [131 show 
that this is equivalent to 

(7) fc<£-l^^VAmax-;^ + 0(l). 

TT 

See Theorem [16] for a more explicit description of eigenvalues with (3 close to some 
/3o>0. 

We can also describe the eigenfunctions in terms of the scattering solutions and 
of the eigenfunctions in £p, with exponentially small errors, see Theorems [211 HH 
and [24l and Corollary [23l for precise statements. 

The restriction Amax < i^i is made to keep the exposition at a reasonable length. 
Our methods are such that a generalization to higher eigenvalues, taking into ac- 
count several thresholds, should be fairly straightforward. 

To identify the numbers bk in Theorem [1] in terms of the data, we recall that a 
quantum graph is given by a metric graph (G with the edge lengths 21^) together 
with a self-adjoint realization of the operator —d^/d^'^ acting on functions defined 
on the disjoint union of the edges, where ^ is a variable along each edge, measuring 
length; such a self-adjoint extension is given by boundary conditions at the vertices. 
We state the boundary conditions obtained in our setting. Denote by 

Ne := ker(Ay^ + ly) 

the eigenspace for v of the Laplacian on the cross section for an edge e. For 
each vertex v let Vy — 0e^„A/'e, then clearly V = 0„ K corresponding to the 
decomposition ([3]). Also, the scattering matrix S{a) is the direct sum of scattering 
matrices Sy{a) : Vy. For each vertex v, Sy{0) is an involution on Vy. Denote 

by Ue the restriction of a function on G to the edge e, and by dnUi,{v) its inward 
normal derivative at the endpoint w of e. 

Theorem 3. The numbers bk in ^ are the eigenvalues of the operator —d^/d^^ 
on the metric graph G with edge lengths 21 e, defined on the space of functions u 
which on each edge e are smooth and take values in J\fe and which at each vertex v 
satisfy the boundary conditions 

(8) (ue(w))e~u G (+1) ^ eigenspace of Sy{0) 

(9) (9„Ue(f ))e~ii G (~1) ~ eigenspace of Sy{0). 
Remarks. 

(1) If, for some edge e, the lowest eigenvalue on is bigger than u, then 
Me — {0}, so Ue has to be identically zero, which means that the edge e 
may be omitted from G. In other words, only the 'thickest' edges contribute 
substantially to the eigenvalue /ife(£), for small e. 
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(2) If Ye is connected then Me is one-dimensional, so it can be identified with 
C. But in general there is no canonical eigenf unction, so this identification 
is not canonical. However, if all Ye are the same connected Riemannian 
manifold (for example, in the case of an e -neighborhood of an embedded 
graph, where they are balls) one may choose the same basis element in each 
Me, SO one may think of u as a complex valued function and Vy ^ Cdcg(D) 
where deg(w) is the degree of v. 

(3) Two special cases deserve to be mentioned: 

• Dirichlet conditions at the vertex v correspond to Sv(0) = —Id. In 
particular, if 5(0) = —Id then the limit quantum graph is completely 
decoupled (Dirichlet boundary value problem on each edge separately, 
without interaction between edges). 

• If all Me are identified with C and Sy{0) has {+l)-eigenspace spanned 
by the vector (1,1,..., 1) then ([8]) says that u is continuous at v and 
([9]) (which is equivalent to {dnUe{v))e^v -L (1, ■ • ■ , 1)^ that the sum of 
normal derivatives of the Ue vanishes at v. This is often called the 
'Kirchhoff boundary condition'. 

(4) Theoremsl^ and\^ show clearly the two ingredients that determine the lead- 
ing behavior of eigenvalues: 

• The determination of S{0) and of the -eigenvalues is a tran- 
scendental problem, depending on the vertex and edge manifolds (for 
example, in the situation of the e -neighborhood of an embedding, the 
angle at which edges meet). 

• Given S{Q), the bk are determined solely by the combinatorics of the 
underlying metric graph. 

In the special case of pure Neumann boundary conditions and of Dirichlet 
conditions with 'small' vertex manifolds the leading behavior of eigenvalues 
is determined by the metric graph alone, see Theorem \25l That is, it is 
independent of the manifolds Xy,Ye. This may be seen as reason why the 
general case is harder to analyze. These cases were known previously, see 
EI], IS] and |2Q]. 

It is known that 'usually' (for example in the embedded situation) there 
are L"^ -eigenvalues t < v. so they can not be neglected. See, for example. 

Using Theorem [3] one can show that for generic geometric data (for example, an 
open dense set of Riemannian metrics on the vertex and edge manifolds) one has 
5(0) = —Id, that is decoupled Dirichlet conditions for the quantum graph at all 
vertices. This will be pursued in a separate paper. 

1.2. Outline of the proof of the Main Theorem. We now give an outline of 
the proof of Theorem [2l where for simplicity we assume that all = 1: Let a: be a 
coordinate on the cylindrical part of X°° and of Xq , measuring length along the 
cylinder axis (going from to A'' from both ends of the cylinder axis). 

The proof consists of two steps: First, we use the spectral data on to con- 
struct approximate eigenfunctions on Xq for large N and conclude the existence of 
eigenvalues and eigenfunctions as claimed. Second, we show that all eigenfunctions 
on Xq are obtained this way. 

For clarity we assume in this outline that all multiplicities (i.e. dimensions of £p 
and Vz) are equal to one and that Tp < v for all p (no embedded eigenvalues). We 
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will call a quantity 'very small' if it is exponentially small as — > oo, i.e. 0(e"'^^) 
for some c > 0. 

Throughout, eigenvalues on Xq will be denoted /z, with jj, = i/ + it ^ > v, 
and spectral values for X°° by A, with X = v + if X > ly. 

Case I: Eigenvalues ^ < v — c: Here, everything is quite straightforward. 

First step: If U is an eigenfunction on with eigenvalue X = Tp < v 
then U decreases exponentially in So if we simply cut off U smoothly 
near x = N, making it zero near x = N, then we get a function on Xq 
which satisfies the eigenfunction equation up to a very small error. A 
simple spectral approximation lemma, Lemma [71 shows that — Aj^-n has 
an eigenvalue very close to A. A spectral gap argument gives a similar 
approximation for the eigenfunction. 

Second step: This procedure can be reversed: If u is an eigenfunction on Xq 
with eigenvalue fj, < — c then it decreases exponentially in x, so cutting 
it off near x — N yields a function on X°° which satisfies the eigenfunction 
equation up to a very small error. Since — Ax°° has purely discrete spec- 
trum near ^, it follows by the spectral approximation lemma that it has 
an eigenvalue very close to /i. Therefore, there can be no eigenfunctions on 
Xq in addition to those constructed in the first step. 

Case II: Eigenvalues G (i^ + e^'^^, Amax]: A basic observation is that any 
eigenfunction u on Xq or generalized eigenfunction on X^ with eigenvalue A G 
{ly, vi) can on the cylindrical part be written as Hu + II^u, where Hit, called the 
leading part, is the first mode in the K-direction. IIu is of the form 

(10) nu = e-*"> + e'"^V, V.i^&V, a■.^^^X^ 

and n^u is exponentially decreasing in x. 

First step: The scattering solutions E = Ea.ip have leading part of the 
form (|10|1 with ^ = S{a)Lp. Since WE is not decaying as a; ^ oo, we should 
only cut ofl[ W^E near a; = TV, in order to obtain a small error term when 
constructing an approximate eigenfunction on Xq from E. Therefore, we 
must require WEa^ip to satisfy matching conditions at a; = that make it a 
smooth function on the cylindrical part of Xq . A short calculation shows 
that this is equivalent to the equation 

(11) M(a, if) = 0, where M{a, ip) := [id - e'°''^^^aS{a)] € V. 

Perturbation theory gives functions Zf and spaces so that the solutions 
of this equation are a = j^Z^{j^) and ip G 'Piijf)- So for these a,ip the 
function obtained from E — E^^^p by cutting of li^E near x = N satisfies 
the eigenvalue equation on Xq with a very small error. Therefore, one gets 
eigenvalues as in ([6|). 

Second step: It remains to show that for large N all eigenvalues fj, d {v + 
e~'^^ , Amax] on Xq are obtained in this way. This is the theoretically most 
demanding part of the proof. Let u be an eigenfunction of —Ax^ with 

eigenvalue /i = v + P'^.P > e~'^^l'^. An argument directly analogous to 
the case [i < v won't work since it only yields that — Ax°° has nonempty 
spectrum near /i, which we know anyway. Rather, we need to show that 
u is very close to some E — E^^ip with (a, ip) satisfying (jlip and with 
a very close to /3. The closeness oi u,E would follow from closeness of 



8 



DANIEL GRIESER 



the leading parts IIu, HE since they control the full function, and this is 
proved in two steps: First, we prove an elliptic estimate which reflects the 
essentials of scattering theory in a compact elliptic problem and implies 
that eigensolutions whose part decays exponentially for x < N (as 11^ u 
does) are close to scattering solutions (whose 11-'- part decays exponentially 
for a; ^ oo) with the same eigenvalue, see Lemmas [18] and [TH This is a 
stable version of the existence and uniqueness of scattering solutions on 
So we obtain some Ep^ipg very close to u. In a second step it is shown 
that (/?, ifo) must be very close to a solution of pT|) : Since u satisfies the 
matching conditions at x = N and u — E is very small, E almost satisfies 
these conditions, so M(/3, ip) is very small; therefore, all that is needed is a 
stable version of the analysis of (llip (which, however, is quite non-trivial). 

Case III: Eigenvalues G [v — c,i^ + e^^^]: While the first step (construction 
of approximate eigenfunction) presents no new difficulties, the second step (proof 
that all eigenfunctions are obtained) is quite delicate. One difficulty is that the 
representation pO|) is not valid at a = (it does not give all scattering solutions). 
There is a straight-forward replacement, however. More serious is showing that any 
eigenvalue € — Cjv) must actually be in (i^ — e~'^^ , v). This again requires a 
delicate stability analysis of the matching condition. 

Special care needs to be taken when the multiplicities are not equal to one, since 
it is not enough just to construct eigenvalues, one also needs sufficiently many. A 
useful tool here is the notion of distance between subspaces of a vector space which 
we recall in Section 

If one is only interested in the case of fixed k as in Theorem [1] then some of the 
proofs can be simplified considerably since the functions remain separated for 
different i then. See, for example. Corollary [231 In particular the proof of Lemma 
[211 simplifies considerably in this case, and Theorem is not needed. 

1.3. Outline of the paper. In Section [2| we introduce the setup precisely and 
some notation. Section [3] introduces the basic analytic tools, most importantly 
separation of variables and distance between subspaces, as well as a basic spectral 
approximation lemma. In Section |4| we recall the facts from scattering theory that 
we need. In Section [5| we analyze which scattering solutions satisfy the matching 
conditions, in particular, equation (jlip . The basic elliptic estimate and some con- 
sequences are proven in Section [51 Theorem [21 and improvements of it are proven 
in Section[71 and Theorem[31in Section[51 where we also discuss some special cases. 

In the Appendix we collect some basic results on one-parameter families of uni- 
tary operators which are needed in the analysis of (fTTj) . Their proofs can be found 
in [7]. 

1.4. Related work. As already mentioned, the Neumann problem was treated in 
m, [6], [H], [H], ig. For Dirichlet boundary conditions, Post jMj derived the first 
two terms of ^ in the case of 'small' vertex neighborhoods, see Theorem [25l In 
the recent preprint [17] Molchanov and Vainberg study the Dirichlet problem and 
show that, in the context of Theorem [1] the ^k{s) — e~'^v converge to eigenvalues 
of the quantum graph described in Theorem [3l this was conjectured in [16], where 
also some results on the scattering theory on non-compact graphs are obtained. 
However, their statements are unclear as to whether the multiplicities coincide; also, 
they do not consider the effect of eigenvalues on X°° or uniform asymptotics for 
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large k. In [T] a related model is considered. The method in the previously cited 
papers is to compare quadratic forms or to show resolvent convergence of some sort, 
and in all cases only the leading asymptotic behavior is obtained. 

Problems of the same basic analytic structure (cylindrical neck stretching to in- 
finity, attached to fixed compact ends; usually with G consisting of one edge only) 
were studied by various authors in the context of global analysis, where they occur 
in a method to prove gluing formulas for spectral invariants. In their study of an- 
alytic torsion (related to the determinant of the Laplacian), Hassell, Mazzeo and 
Melrose [TU] gave a very precise description of the resolvent (in the case of closed 
manifolds, i.e. no boundary, but admitting edge neighborhoods which are not pre- 
cisely cylindrical, just asymptotically), including its full asymptotic behavior as 
e ^ 0, using ideas of R. Melrose's 'b-calculus', a refined version of the pseudodif- 
ferential calculus. More direct approaches were used by Cappell, Lee and Miller [3] 
and by Miiller [18 in the study of the ?7-invariant (for the Dirac operator instead 
of the Laplacian) and by Park and Wojciechowski |l9j. The author and Jerison [8] 
prove a special case of Theorem [2] (where Xq is a plane domain obtained by at- 
taching a long rectangle to a fixed domain, which is required to have width at most 
the width of the rectangle) by a different method (matched asymptotic expansions) 
and use it to prove a result about nodal lines of eigenfunctions; for this, one needs 
to know the asymptotic behavior to second order (i.e. one order more than written 
explicitly in ([2])). 

In the context of this literature, the main purpose of the present paper is to give 
a mathematically rigorous yet straightforward derivation of the limiting problem 
on the graph, allowing various boundary conditions, in a way that admits general- 
ization to similar problems. For example, there are straightforward generalizations 
to higher order operators, systems and Schrodinger operators with potential, as 
long as one has a product type structure along the edges. We use existence of the 
scattering matrix for manifolds with cylindrical ends as a 'black box'. The main 
technical problem is the proof that all eigenfunctions are obtained by the given 
construction; here the neighborhood of the threshold i' requires a special effort 
(this is sometimes referred to as the problem with 'very small eigenvalues' since the 
eigenvalues are exponentially close to ly). 

Notation: As usual constants c, C > may have different values at each oc- 
curence (unless otherwise stated). They depend on the data (the graph, the edge 
lengths, the compact manifolds), but not on N. The scalar product on a space 
X is denoted by {u,v)x — fx ^-i^d the norm by || • 

2. The setup: Combinatorial and geometric data 

The following data are given: 
• Combinatorial data: 

A finite graph G = {V,E), with vertex set V — V{G), edge set E = 
E{G). Loops and multiple edges are allowed. Thus, E may be thought of 
as a multiset of unordered pairs of vertices. If the vertex v is adjacent to 
an edge e, we write w ~ e. A half-edge is a pair (v, e) with w ~ e, and we 
denote by E the set of half edges, except that for a loop e at a vertex v 
the element (w,e) appears twice in E (so formally E is really a multiset). 
Sometimes we denote a half-edge by e if it arises from the edge e. The 
neighborhood of a vertex v is the (multi-)set of half-edges incident to it. 
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It will be useful to think of G as the disjoint union of the vertex neigh- 
borhoods, with the 'ends' of the half-edges glued together appropriately. 
The glueing may be encoded by a map ti, where 

(12) a -.E 

(13) {vi,e)^ {v2,e) 

if the edge e connects Vi and V2 (for a loop e, a maps one copy of {v, e) to 
the other), u is an involution, that is a"^ — Id, and has no fixed points. 

We also assume that a positive number Ig is given for each edge e, to be 
thought of as half the length of e. Denote the shortest half edge length by 

'mill = niin . 



• Geometric data: 

— To each vertex v a compact Riemannian manifold {Xy , gx„ ) with piece- 
wise smooth boundarjo, of dimension n 

— to each edge e a compact Riemannian manifold (Yg , gy^ ) with smooth 
boundary, of dimension n — 1; for a half edge e corresponding to an 
edge e, we set Y- := Ye; 

— to each half-edge (u, e) an isometry h^^e (gluing map) from Y^ to a 
subset of the boundary of X„; we assume that {Xy,gx^) is of product 
type near hv^eiYe), see below; also, for each v the sets /ii,,e(^e) are 
assumed to be disjoint (no overlaps of different edges); 

— partitions of the boundary of each Yf, and of the part of the boundary 
of each Xy which is not in the image of any of these isometries, into 
two pieces denoted by indices D and N (for Dirichlet and Neumann 
boundary conditions) ; we assume sufficient regularity of this decompo- 
sition so that the boundary value problems formulated below (before 
([TB)) ) are well-posed. 

More generally, one may give a pair of non-negative functions a, b (of 
sufficient regularity to make the problem below well-posed) on the 
boundary of each Xy (outside the gluing part) and each Y^, with a,b 
never vanishing simultaneously, to define Robin boundary conditions, 
see below. The D/N decomposition corresponds to a,b being charac- 
teristic functions of a partition of these boundaries into two parts. 

From this data, we define Riemannian manifolds X^, < < oo, with piecewise 
smooth boundary, as follows: First, for each half-edge e = (i;, e) we attach a cylinder 
with cross section Y^, 

Z# := [0,7V) X Y, 

to Xy, using the isometry hy^g. Thus, for each vertex v we get a manifold with 
piecewise smooth boundary 



vN 



where the quotient means that each {0} x 1^ is identified isometrically with a subset 
of the boundary of Xy. 



^for an exact definition of this, see for example for our purposes one may assume that the 
data are such that the spaces defined below are smooth manifolds with boundary. 
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On Z-^ we put the cylindrical Ricmannian metric 

(14) gzlf ^lldx^ +9Y.^- 

By assumption, {Xy,gx^) is of product type near :— h^^eiYe), which means that 
a neighborhood of in is isometric to (0, e) x with the product metric, 
for some e > 0. This ensures that the metrics on and define a smooth 
Riemannian metric on (it also fixes the smooth structure on X^). 
Let 

X^ := U 

(disjoint union) and denote by 

Y:=\Jy^, Z^ ■.= [0,N)xY = \J Z^ 

the cross section resp. the cylindrical part of X^. 

For finite N the pieces X^ are now glued together as prescribed by the graph 
G to give the A/"- neighborhood, Xq , of G. More precisely, a : E ^ E induces 
bijections, also denoted cr, Ij- ^ ^o-(e) (since both of these cross sections are just 
copies of the same Y^), and therefore {N} x 1^ ^ {A^} x ^^(e) for each e G E, and 
then 

(15) X^ := 'W/a 

where X^ is analogous to X^, except that x < N is replaced by a; < A''. In other 
words, Xq is the union of the Xy and cylinders [0, 2N] x Ye for each edge e of 
G, glued together according to the structure of G. (But our x coordinate will run 
between and TV from both ends of the interval.) We also write 

Z^ := W/a. 

Clearly, the D/N decomposition of the boundaries of each Xy and of each Ye 
gives a corresponding decomposition of the boundary of each X^, of Xq and Y, 
and Z^ . The Riemannian metrics define Laplace operators A^n , Aj^-jv , Ay , A^w 
on these spaces, for which we impose Dirichlet boundary conditions on the D part of 
the boundary and Neumann boundary conditions on the N part (but no boundary 
condition at the boundary piece {0} x F of Z^). More generally, one may consider 
Robin boundary conditions, au + bdnU — 0, where the functions a, b on dX^ (etc.) 
are induced by those on the boundaries of Xy, Ye, by making them independent of 
the cylinder coordinate x. 

By d, 

(16) A,.=l-'^ + Ay^. 

Remark. The Riemannian metric (|14p expresses the fact that the cylindrical part 
Z^ has length Nle- Instead, one could use the coordinate ^ — xle, then one would 
get the more standard form 

(17) 9z^ = + 9Y^: ^z^ = -Q^ + ^Y,, £.=xle, 

and < f < leN on Z^ . The (^-coordinate is more convenient for some calcula- 
tions, but in the x-coordinate Z^ is given by the same range of x, < x < N, for 
all e. We switch between these coordinates as is convenient. 
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3. Basics of the analysis 

3.1. Notation. Since Y is compact, the Laplacian Ay has compact resolvent and 
therefore discrete spectrum. Let i/q < I'l < . . . be the eigenvalues of —Ay, with 
finite dimensional eigenspaces Vo,Vi, . . . . We also write 

v := I/O, V Vq. 

Let n be the orthogonal projection to V in L'^{Y), and 11-'- = Id^2(y) — 11. 

Let N G (0,oo]. We decompose any u E Lf^^iX'^) in its 'vertical' 11 and 11-'- 
components, over the cylindrical part : 

u^z" — nu -|- H^u 

where {Ilu){x) :— Il{u{x, •)) and similarly for 11-'-. Here and throughout the paper 
we identify functions on with functions on [0, TV) whose values are functions on 
Y. In particular, 

Uu : [0, N) V, n^u : [0, N) V^. 
Since Y = UegB ^ ^^'^ ^'^ts on each Y^ separately, we have 
(18) L2(y) = 0L2(y_)^ V = ^V^, :=ker(Ay, + z.). 

We write elements ip G L'^{Y) a,s ip = {(p-)-g^^ with G i^(^)- It (fi G V then 
(/3e € 14 Ve. 

The data define two important linear maps L, a : L'^{Y) L'^{y) which restrict 
to linear maps L,<t : V V: 

The map L is diagonal with respect to the splitting (fT8|) and encodes the edge 
lengths 

The map a is defined by the involution a : E ^ E encoding the graph structure: 

(20) l\ "-'^"^j Zv^"^^ }' 
We denote 

(21) := the ± 1 eigenspace of a. 

Since cr is a self-adjoint involution, the decomposition V = V~ is orthogonal. 
For if G V we denote the corresponding decomposition (p ~ + Lp'^ . 
By ([H]) and ((T71), we have 

(22) Az«=i-2__^^^ 

k=0 ^ ' k=Q ^ ^ 

with respect to the decomposition L'^{Y) — ^fc- 



2 '^fc 
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3.2. Matching conditions at a; = A^. We regard functions on Xq as functions 
on X'^ satisfying suitable matcliing conditions at a: = iV. For solutions of tlie 
eigenfunction equation we get: 

Lemma 4. Let u he an eigenfunction of — Ax« with eigenvalue A. 

Then u defines an eigenfunction of — Aj^-n if and only if u extends smoothly to 
X = N and 

(24) ^ 0, {d^u)+ =0, atx^N. 
The upper ± refer to the a -decomposition (|2ip . 

For a > let 

(25) MC„ (u) n^;^^ + ^n{d^u)+^j^ G V. 

Then MCq,(m) = if and only if IIm satisfies the matching conditions at x = iV. 
The particular scaling in (|25|) is motivated by the calculation ([SSjl . 

Proof. This is just the fact that the solution to a second order elliptic partial 
differential equation can be continued across a hypersurface, as a solution, if and 
only if at the hypersurface it is continuous and its normal derivatives match. Set 
$ = u^=N, * = id^u)^^N- Then = i(4> - cr$) and = i(* + cr^'), so ((24)) 
means $e = *i'<j(e)i — ~^a{e) for all e, that is, continuity of u and of d^u at 
x = N in Xg. □ 

For a function on X^ , when we write Ux=n we always assume that u extends 
smoothly to a; = A^. 

3.3. Separation of variables. The following simple lemma is basic to all the 
analysis. 

Lemma 5. LetQ<N < oo. Let u e C^iZ^) satisfy (A^iv + \)u = 0. 

(a) The leading part of u has the form 

(26) IVu = e-^^^if + e^^^^V if >^ < 

(27) Hm = e~*"V + e'"^V' if\>v, \ = v + a^ 

(28) Ilu = ^ + ii> «/ A = I/, 

wif/i t/" G (Replace ^ by xL to express this in terms of x.) 

(b) Lf X<vi then 



oo 



(29) n^u = e-^^^^k + e^^^^Pk, Vk,i^k & 14 

fc=i 

(c) ip,ip in (a) and ipk,tpk in (b) are uniquely determined by u. Lf N — oo and 
u is polynomially bounded as ^ oo then tp = in (|26p and V'fc = for all 
k in (pg)) . 

(d) Assume N < oo and u extends to a solution on Zq , or N = oo and u is 
polynomially bounded. Then there is a constant C such that for all M < N 

we have 

(30) Wu^^mWy + \mu)^=M\\Y < Ce-'''\\u\\zi 
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if X < ly, where c = \/ v — X /min> 
(31) \\^^u^=m\\y + ||(n^a^^),=M||y < Ce-=*^||n^ii||2i 



if X < i^i , where c = y/i/i — X Imin ■ 
(e) If N < oo and u is an eigenfunction on Zq then, for the representation in 

(EH), 

oo 

(32) II Y^ivkf/^V^^^MW < Ce-2^^||n^u,=o||Y 



if X<vi for any c < ^/vi - Xlmin- 

Proof. For fc G No let 11^ : L'^{Y) Vk be the projection and let Uk = IlfcU : 
[0, A^) Vk- By Ufc satisfies the differential equation 

(33) (^ + A - i.fc)wfe = 0, 

and this has the solutions given by the formulas in (a) and by the summands in (b). 
Clearly, Uk and therefore (p, ip, ipk,''pk are determined by u, and if u is polynomially 
bounded then so is Uk = Ilfctt for each fc, so (c) follows. This immediately gives 
(d) in the case N = oo. If N < oo and u extends to a solution on Zq then u may 
be regarded as function [0,27V] L'^iY). We now use the coordinate x S [0, 2iV]. 
If i^fc > A then one can express the solution of ((33|) by its values at a,b, for any 
< a < 6 < 2A^: Write Pk = Vi^k - XL, then 

sinh(.T - a)Pk . , sinh(6 - x)l3k , 

smh(o — ajp/c smh(o— a)pfe 

For any a € (0, 1) and h G (2A^ - 1, 2N) one obtains easily at x = A/ < iV: 



lkfe(M)||v, < Ce-^^'"""*^ (hfe(a)l|y. + ll"fe(6)|kJ • 

Summing over fc = 0, 1, 2, . . . in the case A < yields the estimate on Ux=m in 
(150)) . The estimate on the derivative and estimate (|3ip are obtained similarly. 



(e) fo llows similarly to (d), using II-'-m^j^o = J2k Vk+i'k, 'n^Ux=2N = Efc e ^"'^ ^^(/?fc- 
gVi^fc-AAf^^ and U2,=27V = crw^^o- □ 

It is clear from the proof that the behavior of solutions is different for A > i^i . 
For the purpose of this paper, we will always consider X < vi. Define Amax,Q!max 

by 

(34) Amax = ly + a^ax7 ^max £ (i^, vi) arbitrary, 

then the estimates will always be uniform for A < Amax- 
We define the 'boundary data' of a function u on Z^ as 

(35) w° nu|2,^o, ""^ n(9^u)|j,^o- 

Clearly, if u is an eigenfunction then {vP,u^) determine IIu uniquely: Instead of 
the representation ([26]) . ([27|) . ([28]) for the leading part of an eigenfunction we use 
the basis 

fcosVA^C (A>z.) f'^^^^ (^>^) 

[coshV^T^e (A<i^) I ''°vS^^ (A<i^) 
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of solutions to {d^ + A — iy)v = 0. This is more useful than the exponential basis 
for A near i' since it depends analytically on A even at A = i^. For a solution u of 
(A^jv + X)u = we have 

(36) Ilu = Cxu'^ + Sxu\ 

3.4. Distance betw^een subspaces and spectral approximation. It will be 
convenient to use the notion of distance between subspaces. Although it is standard, 
we recall its definition and basic properties here. 

If V, W are closed subspaces of a Hilbert space H then define 

dist{V,W):^ inf '^'ffl 
uev\{o} \\u\\ 

where dist(u, M^) — inf{||u— w\\ : w £ W}. dist is not a distance function since it 
is asymmetric (it does satisfy the triangle inequality, however). 

It is elementary to check the following properties (see [12], for example): 

Lemma 6. a) dist(y, W) = \\PwPv ~ Pv\\ where Py, Pw a'^e the orthogonal pro- 
jections to y, W , respectively. 

b) dist(V^, W) < 1, and < 1 ijjllw\v -V is injective. 

Here, Tlw\v the restriction of Pw to V. In particular dimV^ < dimVF if 
dist{V,W) < 1. 

c) If dist{V,W) < 1 and dist{W,V) < 1 then Pv\w ■ V W is an isomorphism 
and dist(y, W) ~ dist(T4^, V). In this case we write 

distsymm(^, W) dist(F, W) = dist(M/, V). 

d) IfVi,i^ I, are pairwise orthogonal then 

(37) dist(0 V,, W)<Y, dist(^^^, W). 

i i 

e) If Ui, i = 1,...^K are pairwise orthogonal and dist(span{ui}, T4^) < (1 + 
dimVF)"^ for each i then K < dimM^. 

(Proof of (e), for example: Let K' — dim W. UK > K' then dist(span{Mi, . . . , uk'+i\^ W) < 
(K' + 1)1 {K' + 1) = 1, so K' + l< K' by (d) and(b), a contradiction.) 

We will use the following standard spectral approximation lemma. We include 
a proof for completeness. 

Lemma 7 (Spectral Approximation Lemma). Let A he a selfadjoint operator in a 
Hilbert space TL. For an interval I G R let Eigj(A) be the spectral subspace of A 
corresponding to the spectral interval I . 

Let Aq G K, e, (5 > 0. IfWd Dom(74) is a linear subspace satisfying 

(38) ||(A- Ao)m|| < £||u|| yueW 
then 

(39) dist{W,Eig^Xo-SM+sM))<'^- 

In particular, if e < 5 then dm\'EAgi^-y^^_g -^^^g-^{A) > dimW^. 

An important consequence is that existence of spectral gaps implies good ap- 
proximation of eigenfunctions: If A has no spectrum in {A : e < |A — Ao| < <5} then 
W is e/(5-close to Eig(-;^jj_g (A). 
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Proof. We may assmiie Aq = 0. Let A = j \ dE\ be the spectral resolution of A. 
Let Ps be the the projection to Eig(-_^ If m e Dom(A) is arbitrary then 

(40) 5''\\u~ Psu\\^ <\\Auf 

since J|^l^^- d(w, £'a(m)) < ^^^^y^\^d{u,E\{u)) < j^\^d{u,E\{u)). Therefore, if 
w eW then \\w — Psw\\ < ^\\w\\, which was to be shown. The last claim follows 
from Lemma |B)d). □ 

4. Scattering theory background 

Here we collect the facts from scattering theory that we need. A good reference 
for this material is [9j (this may not be the earliest source). 

Scattering theory is the spectral theory of an elliptic operator on a non-compact 
space which has a 'simple' structure at infinity, that is, is asymptotically equal to an 
operator on a space for which the spectral decomposition may be written down fairly 
explicitly. Among the central goals of scattering theory are the determination of 
the absolutely continuous spectrum and of generalized eigenfunctions correponding 
to it. 

We will need fairly explicit spectral information on the operator — Ax°°- Since 
X°° is cylindrical at infinity (that is, outside a compact subset), we take as explicit 
model — A^oo, where the prime means that we impose Neumann boundary condi- 
tions at a; = (Dirichlet would be equally possible), in addition to the boundary 
conditions coming from the D/N decomposition (resp. Robin data) of dY, in order 
to make the operator essentially self-adjoint. The spectral theory of —A'^^ is easy 
to obtain, using separation of variables, i.e. ((23|) . We use the ^-coordinate on [0, oo) 
for simplicity. It will be replaced by xL at the end. 

,2 

Since —-j^ on [0,cxd), with Neumann condition at ^ = 0, has absolutely contin- 
uous spectrum [0,cx)), with generalized eigenfunctions cosa^ corresponding to the 
spectral parameter a^, a > 0, the decomposition (|23|) shows that the absolutely 
continuous spectrum of the model is 

speCabs(-A2^) = IJ K,oo) = K,oo) = hoo), 

A;>0 

and that the functions 

cos{\/ A — i^k for k,ip satisfying J^fc < A, ip € Vk 

span the generalized eigenfunctions of — A'^oo with eigenvalue A. That is, each such 
function U satisfies —A'^oaU = XU and is polynomially bounded as ^ — > oo, and 
any function with these properties is in the linear span of these functions. 

The scattering theory for spaces with cylindrical ends shows that this picture 
carries over to — Ajf°o, except for the possible appearance of discrete spectrum, and 
a 'phase shift' and exponentially decaying error term in the generalized eigenfunc- 
tions. 

First, we have the description of the spectrum: 

Theorem 8. ([9\) The operator —Ax'=° is essentially self-adjoint on C^{X^). Its 
unique self-adjoint extension (still denoted — Ax°° ) has the following properties: 

a) The pure point spectrum is a discrete subset of [0, oo). 

b) The singularly continuous spectrum is empty. 

c) The absolutely continuous spectrum is [j/, oo). 
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Discreteness includes finite multiplicity. 

We now discuss the generalized eigenfunctions (or 'scattering solutions') E, 

(41) -Ax-^E^XE, \ = v + a^ 

for spectral values A € [v,vi). One could also consider values v > vi, this would 
mean using higher scattering matrices. 

By Lemma El E\z-=^ ^ HE + li^E where E has the explicit form (gT]) or (QSl) 
and n^E is a sum of exponentially decreasing and exponentially increasing terms. 
By polynomial boundedness, the latter must vanish. 

Theorem 9 ([9]). 

a) (\ > v) Assume that A G is not an L'^ -eigenvalue of —Ax°^ ■ For each 
If d V there is a unique bounded solution E = Ea^ip of (I4ip satisfying (with 
a — \/\ — v) 

(42) ni;„,^ = e"*"V + e'"V ^or some ip' G V. 

b) (Scattering matrix) This defines a linear map 

S{a):V^V, ip^~^ip', 

called the scattering matrix. S{a) G End(Vb) is holomorphic in a and extends 
meromorphically toaGC\{Q!GM; |a|> — v} , holomorphic for real a. 
Furthermore: 

i) S{a) is unitary for a real. 

ii) S{a)S(—a) — I for all a. 

ill) S{Q) is an involution, and S"(0) commutes with S{Q). 

c) (\ — v) The space Sa — {Ea,tp ■ V? G V} depends holomorphically on a for a as 
in a) and extends holomorphically as in b). For a — Q it is described as follows: 

Let V± be the ±1 eigenspaces of S{0). For all G V+, G V- there is a 
unique generalized eigenfunction E — -Eo.'i>+.*_ G satisfying 

(43) n£;o,<i>+,*_ = $ + $ = *+ + ^5"(0)^'_, 

Here, ^^(O)*- G y_ . 
The functions in (a) and (c) are all the generalized eigenfunctions, up to addition 
of possible Lp' -eigenfunctions. 

Note that if v -\- is an L^-eigenvalue, the function Ea,i^ resp. -Bo,$+,*_ is 
not uniquely determined by fixing its leading part, since the leading part of an 
L^-eigenfunction is zero. However, it is determined by the additional requirements 
that it lie in Sa and that this space depends continuously on a (for a 7^ this 
means that a i-^- E^^ip depends continuously on a). 

Proof. This is mostly standard. Note that b)iii) follows from ii) by setting a = 0, 
and from differentiating ii) at a = 0, which gives S"(0)S'(0) - S'(0)5'(0) = 0. This 
implies that 5'(0) preserves V±, and therefore the last claim in c). Equation (l43l) 
will be explained below. □ 

In our context, we use the variable x, where ^ = xL. Therefore 

(44) IVEa,^ = e-'°''-'^ip + e'"^-^S'(a)(/7 

(45) n£;o,<i>+,^._ = ($+ + ^5'(0)*_) + xL^-, 
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respectively. 

Theorera [9] can be reformulated in terms of 'scattering subspaces'; we will use 
this since it allows a more uniform treatment of the cases a ^ and a = 0: 

Recall the notation w'^ = Ilux=o, = ^d^u.x=o- Define the 'scattering subspace' 

for a e [-Qfinax, Qfmax] by 

(46) Ca {{E°,E'^) : (A^- = 0, E polynomially bounded} cVxV. 

Since HE determines and is determined by E^), part a) of Theorem [9] may be 

reformulated as follows: 

(47) 

-Cq = {( (/ + S{a))ip, —ia{I — S{a))ip ) : £V} for a 7^ not an L^— eigenvalue, 
and to each (i?°, E^) ^ Ca there is a unique scattering solution E with eigenvalue 

In particular, Ca is a dim ^-dimensional subspace oiV y.V for these a; since 
Ca = IIiSq it extends analytically to all |a| < ^Jvi — v. However, in general the 
continuation to a = cannot be obtained by setting a = on the right of ([17|) 
since this yields V+ x {0}, a space of lower dimension (unless 5(0) — Id), so this 
cannot be Lq. In other words, the parametrization of £q by does not extend 
uniformly to a = 0. Therefore, we will also need another parametrization which 
works (and is analytic) also at a = 0: 

Lemma 10. For a > 0, 

= {(£;", E^) ■.E'^^p+ + f 5'(0)p_ + aR°{a)p 

(48) = p_ - iaS"(0)p_ + a^R\a)p, 

for certain families R'^{a),R^{a) of endomorphisms of V , depending analytically 
on a for \a\ < amax- 

If S{a) = S'(O) for all a then Ca = {{p+-P-) ■ P G V} for all a. 

Proof Write S{a) = 5(0) + aT + a^R{a), T 5'(0), by Taylor's formula. By 
(|Tfl). {E^, E^)cV xV is in (for a > 0) iff, for some ip(^V, 

E" = {I + S{a))ip = 2lp+ + aTif^ + aTip+ + a^R{a)ip 

E^ — ~ia{I — S{a))ip — ^2ia(p^ + ia^Tip^ + ia^T^p-^- + ia^ R{a)<p. 

Here we used 5(0)(y5 = (^+ — and [T{p)± = Tip±. Write 

(49) p+ = 2ip+, p_ = -2iaip-, 

then this becomes (|48p , with suitable R^,R^. The last statement is clear from this 
derivation. □ 

In particular, 

(50) /:o = {(p+ + ^5'(0)p_,p_) : peF}, 
and this explains (|43l) . 

5. Scattering solutions and matching conditions 

In this section we analyze for which scattering solutions E the leading part HE 
satisfies the matching condition at x = iV. 
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5.1. The case \ = v. 

Lemma 11. Let $+ e V+, G VI. If N > 27^||5'(0)|| then n£;o,<i>+,*_ , given 
by (|45p . satisfies the matching conditions (|24p iff 



(51) CT$+ = $+, *_=0. 
Proof. The matching conditions ([24| for ni?o, are 

(52) (/-ct)($+ + ^S"(0)«'_ + ArL*_) = 0, (/ + o-)*_=0. 

This is clearly satisfied if ^ holds. Conversely, if ([52]) holds then $+ + |S"(0)^'_ + 
NL'i>- and are orthogonal, so 

(53) = ($+ + ^5'(0)*_ + 7Vi4»_^ = (ls"(0)«'_, + Ar(i*_, ^f^), 
and Cauchy-Schwarz gives 

||^5'(0)M'_|l.||vl/_|| <l||5'(0)|l.||vl/_||2. 

Since Unll^-lP < (L*-,*-), this and dSS]) imply = if iV > 2i^||S"(0)||, 
and then ([52]) gives a$+ = □ 

5.2. The case X > v. First, we express the matching conditions in terms of the 
scattering matrix. 



Lemma 12. Let a > 0, if £ V . Then Ili^ct^,^, given by (|44p . satisfies the matching 
conditions (|24p 

(54) {I-e"'^^^aS{a))ip = 0. 
Proof Let V = S{a)(p. Then by ([21]) 

a 

(55) ^ g-»aiVL ^ e^'^^^^A]- + [(^ - e2'"^^^]+) 

= e-''"^^((p-e2*"^^a^). 

Here we used that aL = La. (This reflects the fact that two half-edges corre- 
sponding to the same edge have the same length.) Ili^Q.^^ satisfies the matching 
conditions iff MCa^Ea.^) = 0, so the claim follows. □ 

We now analyze the solutions {a,ip) of equation (j54p . and in particular their 
asymptotic behavior as — s- oo. It is convenient to introduce the rescaled variable 

z := aN. 

Let 

U{z,a)=e''^^(jS{a) 
for z > 0, where a € [0, amax]- We need to study the zero set 

(56) Z = {(z, a) e (0, c3o) x [0, a^,^] : det(/ - U{z, a)) = 0} 

and to each (z, a) d Z the eigenspace ker(/— C/(z, a)). Intersecting Z with the line 
z = aN then gives the solutions of ((54|l . 
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The U{z, a) are unitary operators on V. Essential for the sequel is the 'mono- 
tonicity' in z: 

(57) -^U-' - 2i > 0. 

I oz 

We begin with the case a = 0, which should be regarded as the limit N ^ oo. 

Proposition 13. The set of z for which {z,0) G Z consists of a sequence < zi < 
Z2 < ■ ■ ■ ^ OO. IfVzi denotes the kernel of I — U{zi, 0) then 

(58) y diniT'^. - — A + 0(1), A ^ oo. 

i:Zi<A 

This is obvious in the case L — II, i.e. if all edges of the graph have equal lengths, 
since then Vz is just the eigenspace of (jS{0) with eigenvalue e~*^^'. 

Proof. By monotonicity, equation (|57p , we can apply Lemma [26] to the unitary 
family z i— > C/(z, 0) with D{z) — 2L, and this gives the result. □ 

The structure of Z and the eigenspaces is then given as follows. 

Proposition 14. Let zi, Z2, . . . be as in the previous proposition. For each i € N 
there are € N and pairwise different real analytic functions zf , p = l,...,ri, 
defined on [0, amax], such that 

(59) z?(0) = Zi for each i,p 
and 

(60) Z = |J{(zf (a), a) : a e [0, a^,^]}. 

i,p 

There is a constant Co such that 

(61) |^2f(a)|<Co for alii, p, a. 

Furthermore, for each i, p there is a real analytic family of orthogonal projections 
Pi (a) onV,aG [0, amax]? such that for each (z, a) G Z 

(62) P[{a) — the projection to ker(/ — [/(z, a)). 

Thus, Z is the union of graphs of functions of a with bounded derivatives, see 
Figure [21 The 'non-linear eigenvalues' z^ at a = may bifurcate into various z^ as 
a increases. The sum in (j62[) is over the various branches that meet at (z, a), so for 
almost all (z,a) it has only one term (in particular P^ is uniquely determined). 

Proof. Theorem [29] shows that Z is everywhere locally a union of graphs of analytic 
functions, and these can be patched to functions on all of [0, amax]- The theorem 
gives also the projections. Applying the theorem at (zi,0) gives 

It remains to prove (|6ip. If z{a) = zf(a) and 1^9(0;) G RanP/'(a) is normal- 
ized and chosen analytic in a then differentiating U{z{a),a)Lp{a) = '^{a) gives 
z'{dzU)ip + {daU)(p + Uip' — (fi' . Taking the scalar product with (p one obtains, 
since {U(p',(p) — {(p',U^^(p) ~ {ip',Lp), that z' {{dzU)ip,ip) + {{daU)(p,ip) — and 
hence, using ip ~ U^^ip and (|57)) . 

(e»^2^|c75'(a)y>,y) 
(2L^,^) 
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Figure 2. The solutions of det(/ - U{z, a)) = 



which is uniformly bounded as claimed. □ 



We now study the solutions of (|54|) and for this the intersection of Z with the 
line z — aN , for fixed N. The following is quite obvious from Figured 

Theorem 15. (a) Fix arbitrary i,p. Then there is a unique point 

(63) i'^i -^i n) ~ intersection point of z = (a) and z — aN 



whenever N > No^i Cq + with Cq from (j61|) . 

Furthermore, is real analytic in , i.e. zf^ — Z^{jj) for a function 
Zl which is real analytic on [0, {No,i)~^]. Also, Z^{0) = Zi. 
(b) The pairs {a,ip), < a < amax, for which IlEa^,p satisfies the matching 
conditions l\2A\i are given by 

(64) a = ^ ^z^^ 

(65) y^eV^^N-^ RanPfia) 

for i — 1,2,... satisfying N > A^o.i o,nd p = 1, . . . , r^. The P[ are defined 
in Proposition \14\ 

Proof, (a) Fix i,p. The z-coordinates of the intersection points (|63p are the so- 
lutions of z = z'^{tz), where t = j^, that is, the zeroes of the function d{z) = 
z - z^{tz). We have d{0) = -z, < and d{^) = ^ - (ofmax) > for 
t < imax (^^o,i)"^ since zf (a,„ax) < Zi + Coa^ax by integration of the bound 
(|6ip . Therefore, d has a zero z for each t < fmax- The zero z is unique since 
d'{z) = 1 - t{z^y{tz) > for t < 1/Co (which is satisfied for t < t,„ax)- Clearly, 
for t = the solution is z = by ((59|) . and the inverse function theorem gives the 
analytic dependence on t. 

(b) By Lemma [T2| IlEa,ip satisfies the matching conditions iff ([54|) is satisfied. 
The claim then follows from Proposition [Ml □ 
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A more precise analysis yields the uniform behavior of the functions with 
respect to i — > oo: Z^(t) = Zi + h'^{tzi) for functions ft,^ which vanish at zero and 
have bounds on their derivatives independent of i. We will omit the proof since 
we don't need it. Instead, we prove the following slightly weaker consequence of 
Theorem [15] about special regimes of solutions: 

Theorem 16. Given aa S [0, Ofmax] cmd S > 0, C > 0, the values of a in (j64p in 

the interval \a — ao\ < CN^^ satisfy 

where — z'^{aQ), that is, b — b^ are the solutions of 
(67) det(/-e'''2-E'cr5(ao)) = 0, 

and only the range |6f — A^aol < CN^~^ is considered. The constant implied in 



66|l is uniform in i, p, N and S. 



(68) "-l^ + 0(^) 



In particular, for ao = 0, (5 = 1 we get: The solutions a with a < C/N are of 
the form 

f + < 

This is already clear from Theorem [TSb) . 

Proof Let t=l/N. If a = (t) satisfies \a - ao\ < Ct^ then \z1{a) - (ao)| < 
CoCt^ by (mi). From a = tz{{a) it then follows that \a ~ tb1\ < CoCi''+\ that is, 
(|66l) . Finally, this estimate and \a - ad < Ct^ imply \tb1 - ao\ < Ct^{l + Cot), 
and this yields the restricted range for 6f . □ 

We will also need the following stable version of Theorem [TSb). 

Theorem 17. Let N > 0. Assume a e (0,Q;max) md (p d V are such that HEa^^ 
satisfies the matching conditions up to an error 

(69) \\MCo.{E^,^)\\v <6Mv, 
for some 6 > 0. 

a) Then there is i, p with 

(70) |a-<^|<-^A 

b) Furthermore, there is C > such that, if S < C^^ and N > C, then the sum 
'Pa,s,N = ®pT''p,N, with (3 ranging over \a — (3\ < N^^\/S, is direct and, if P is 
the projection to Va,s,N , then 

(71) y-P^Wv < ,5i/(2+2dimy) ii^ii^ 

Proof By ([55]) we have e'"^^ MCc{Ec,,p) = (/ - e'"^^^aSia)) ip, so §9^ implies 
- e'°''^^^aS{a))ip\\ < S\\ip\\. Write z = aN and 

(72) u{z)^e^^'^aS{^), 
then this reads 

(73) \\{I-Uiz)M<6M. 
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Now z U{z) is a monotone unitary family in the sense of ()130|) . (I13ip since, with 
a prime denoting derivative in z, 

-U'U-^ = 2L+ J-US-^S'U-^ 
I iN 

where U is evaluated at z and S at z/iV, and for N sufficiently large and a = z/N 
bounded this is positive with bounds as in (|f 3f p . where dmin can be taken as Zmin- 
Therefore, we can apply Lemma E51 with e — 5. ([70)1 follows from (|135p (remem- 
ber a = z/N) and ^ from (fT36l) . □ 

6. Elliptic estimates 
We first state some fairly standard elliptic estimates on the 'compact part' 

X° U X, 

of X°° . Since we consider the scattering theory results in Section S] as a black 
box in this article, we derive them from those results. In a more thorough and 
systematic treatment, they could be derived directly from the theory of elliptic 
boundary value problems and then used to derive the scattering theory results. 
However, the boundary value problem is slightly non-standard since it involves a 
non-local (pseudo-differential) boundary operator. 

The boundary dX^ splits in two parts: The part where the cylinders are 
attached, which we may identify with Y (which we sometimes also call {x = 0}), 
and the complement dX^ \ Y. At the latter, we have boundary conditions given by 
the D/N decomposition (resp. Robin data). At Y we will now impose boundary 
conditions motivated from the scattering theory. 

For u e C°°(X") and X < i^i let 

B^{X)u := n^d^u^Y - Q(A)(n^W|y), 
where for (^^ G 14, fc e N, 

C30 OO 

k=l k=l 

if the sums converge. Q{X) is the n-*^ part of the Dirichlet to Neumann operator for 
exponentially decreasing solutions of (A^oo -f- X)u = 0, see ([29| . Thus, B±{X)u = 
iff n^u has no exponentially increasing part, where u is the unique function on Z°° 
satisfying (A^oo + X)u = and having the same value and normal derivative at 
X = as u at y. 

Consider the operator A^o with domain Dom(Axo) C H'^{X^) defined by the 
D/N boundary conditions (resp. Robin data) at dX'^ \ Y. In order to obtain a 
selfadjoint extension of A^o we need in addition to impose boundary conditions 
at Y. In addition to the condition B±{X)u = we need a condition involving 
M° := nu|y, :— Hd^u^Y- It is well-known and easy to check that selfadjoint 
boundary conditions correspond to subspaces C G V xV which are Lagrangian, i.e. 
such that dim C = dim V and 

(m°, v}), (w", i;i) e £ =^ (li", f 1) - (w", u^) = 0. 

Thus, for C Lagrangian and any A the operator A^o is selfadjoint on the domain 
{u G Dom(Axo : B^{X) = 0, {vP,v}) G £}. 
We have the following elliptic estimates. 
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Lemma 18. Let a G [0,Q;max]- Assume X = v + is not an eigenvalue of 
—Ax°° ■ Let Ca he the scattering subspace ()46|) . and let C'^ be a Lagrangian subspace 
ofVxV such that n £^ = {0}. Also, let Pc^.c -VxV^Ca be the projection 
along CJ^. There is a constant C so that if u G C°°{X'^) satisfies the D/N (resp. 
Robin) boundary conditions at dX^ \ Y and 

(A + A)u = / m X° 

(74) B^{\)u^g atY 

Pc^,C'M'y)^h 

then 

(75) ||u||j^2(jf0) < C( + 11511^1/2(1-) + \\h\\vy,v )• 

// a is varied and the CJ^ depend continuously on a then the constant C can be 
chosen independent of a. 

See Section 17.41 for a replacement in case ly + is an L^-eigenvalue. 

Proof. If u is a solution of the homogeneous problem, i.e. f ^ 0, g ^ 0, h ~ 0, then 
g = implies that u extends to a solution of (A + X)u = on X°^' with bounded 
n-*- part and hence that u is a scattering solution. Then (m*^, u^) G Ca by definition 
of Ca- h = implies {u°,u^) £ C'^, hence {u°,u^) = since Ca r\ C'^ — {0} and 
therefore u = by uniqueness of scattering solutions. 

Therefore, the map u ^ ((A + A)u, -B_l(A)m, P£„,£^ w^)), Dom(Axo) 
L^(X°) X H^/^{Y) X {V xV) is injective. It is also surjective, since ior g = 0, h — 
one has a solution u for any / by the Fredholm alternative (since the operator A^o + 
A is selfadjoint on the domain which consists of those u satisfying homogeneous 
boundary conditions, has closed range and is injective), and arbitrary g,h can be 
removed by replacing m by m — u, where v is any function in Dom(Axo) satisfying 
B±{X)v = g, Pc„,£^(v°, w^) — h. V exists by standard arguments, since X^ is 
of product type near Y by assumption. Since Dom(Axo) is complete with the 
i?^-norm the open mapping theorem gives (j75[) . 

To show that C can be chosen independent of a it suffices to show that it can 
be chosen locally uniformly with respect to a. This can be proved as follows: Fix 
ao and let Co be the constant for a = ao. Suppose u satisfies (1741) for some a near 
OfQ. This can be rewritten (A + Xo)u = / + (Aq — X)u, B±{Xo)u = g + {B±{Xo) — 
Bj^{X))u,Pc^^^C'^^iu'',u^) = h+{Pc^^^,c'^^ - P£„,£j(u",ui). Estimate ^ with 
these data yields ||u||//2(xo) < Cq{ ||/||l2(xo) + \\g\\m/^(Y) + £\\u\\m(xo)), where 
e ^ as a ^ ao (independently of u), since all operators on the right in 
depend continuously on a. For e < l/2Co the last term can be absorbed into the 
left hand side, and the claim follows. □ 

Using ((75)l with u = E a. scattering solution we get 

(76) \\E\\h^(^xo) < C'\\{E° , E^)\\vxv for any scattering solution i?. 

When applying Lemma [15] we will need the following estimate which shows that 
the exponentially increasing part, B±{X)u, of an eigenfunction on Xq is very small. 

Lemma 19. Let (Ajj-n + /i)M = 0, fi < i/i. Then 

(77) \\BAl^Mm/^Y) < Ce--^\\u\\H^xo). 
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Proof. If U-^u = J2kLi e-^'^^-'^^ipk + e^'^'=~^«Vfc with (pk,iJk & Vk then B^{ii)u = 
Sfc°=i V^fc ~ mV"*;- Since Ay is an eUiptic operator of order 2, with eigenvalues v^, 
it foUows that \\B^{^)u\\hi/2(y) < II Sfcli(^fe)^^''V^fc ~ f^i'kWL^iY), and then ((321) 
together with the trace theorem gives the claim. □ 

We will need that scattering solutions for different spectral values whose leading 
parts satisfy the matching conditions are almost orthogonal on X^: 

Lemma 20. Let a, a' £ An and ip e Va,N,'fi' G Va'.N- If ci ^ a' then the 
restrictions of E ^ Ea^ip,E' = Ea\ip' to are almost orthogonal, i.e. 

\{E,E')xN\<Ce--^\\E\\xA\E'\\xr.. 
Proof. With \ = V + , \' = v + (a')^ we have by Green's formula 
(A — A')(-£/, E'^x^ ^ (-^-^5 I^')x^ — (-^7 X' E')x^ ^ — (A-£/, E')x^ ^ (-^5 ^E') x^ 
= {Ex^x,d(^E'^^j^)Y - (d^Ex=N,E'^^!^)Y 

Since HE, HE' satisfy the matching conditions at a; = iV, the 11 part of the lat- 
ter scalar products vanishes, and only the 11^ part remains. Writing U^E = 
e'^^^^^^^^ipk, tpk €Vk, and similarly for 11-'- i?', we get 



(A ~ X'){E, E')xN - ^(v/^^^ - v/^;:^)(e-^^^^v'fc, e-^^^^^y.;,) 

A.— 1 

00 

= (A' - A) v^;;^)-i(e-^^^^^,,e-^^^ 

k=l 

and since the latter sum is bounded by 




Ce-^"^ JE ll^fcll' JE ll^^ll' = ^^"''^lin^^-oll^lin^^-oll^' 



and since ||i?2;=o|| < C||i?||xo by elliptic regularity, the claim follows. □ 



cases: Eigenvalues on Xq arising from L^— eigenvalues on X°° below the essential 
spectrum; eigenvalues on X^ arising from the continuous spectrum on X°^ but 



7. Proof of the Main Theorem 
In this section we prove the main theorems. We treat separately the following 

N 
G 

away from the threshold, i.e. bigger than + e"'^^ for suitable c > 0; eigenvalues 
on X^ arising from the threshold. 

In each case, the proof proceeds in two steps: In Step 1, we construct the eigen- 
values on Xq from approximate eigenfunctions constructed from the (generalized) 
eigenfunctions on X°°. In Step 2 we show that in this way all eigenfunctions are 
obtained. 

We assume at first there —Ax°° has no L^-eigenvalues in [ly, Amax]- The modifi- 
cations needed in case there are such eigenvalues are described in Section 17.41 

We use the following notation: For a selfadjoint operator A and / C K let 
Eigj(A) be the spectral subspace for the spectral interval /. If A has only discrete 
spectrum in /, this is the span of the eigenfunctions of A with eigenvalues in /. 
Also, let 

Eig/,jv :=Eigj(-Ax^). 
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We will construct approximate spectral subspaces Eig^''^, defined in each case sep- 
arately, using a cutoff function defined as follows. Choose % € C°°(M) satisfying 
= 1 for a; < — i and x{^) = for x > — ^ and set 

(78) XN{x)^xix-N). 

Thus, if u is a function on then xnu equals u for x < N — 1 and is identically 
zero ior X > N — ^. 

Note that one could not hope to approximate individual eigenfunctions on Xq 
by approximate eigenfunctions if the eigenvalues lie very close together. 

7.1. Eigenvalues arising from eigenvalues on X°°. For / C (0, j^) let 

(79) Eigf P span{xiv" : u e Eig,(-Ax-)}. 

Theorem 21. For any cq > there is Nq > such that for N > Nq the eigenvalues 
jJL < V — of — Aj5f jv lie within Ce~^^ of the -eigenvalues of —/S.x°° ■ 
For each L'^ -eigenvalue \<v o/— Ax°° we have 

(80) dist„(Eig, Eig^A},iv) < ^e-^^, / = (A - Ce^^ , A + Ce-^^). 



Here c = \/v — A^ax where A^ax «s the largest eigenvalue o/— Ax°o less than v. 
Proof. Step 1: Show that the approximate eigenfunctions are actually such: 

yapp 
'{A}, A" 



(81) dist(Eig^PP ,Eig,^)<l. 



In particular, — Ajs^-w has at least dimEig^^^? many eigenvalues in /. 



Proof: Let W = EigJ^^ For w = xnu € W one has 



={A},Af 

Aw xa(Am) + 2VxArVw + (AxAr)w = -AxATW + 2L-'^x'NdxU + L-^XnU, 

and since Xat; Xn supported in {N — 1, N), one obtains from pO|) . applied with 
M G {N-1,N), that 

(82) ||(A;,«+A)HIx« <Ce-^'^||«;||x«. 

Since W C Dom(Ajfjv) we may apply the Spectral Approximation Lemma[7]to the 
operator A — — Aj^-n, with Aq = A and e — Ce^"^^, S — 2e, and this gives ([5T|) . 

Step 2: Show that any eigenvalue /i of — Aj^-n is in some / = (A — Ce'^^, A+Ce^^^) 
and that 

(83) dist(Eig,,^,Eig^PP^)<Ce-^^. 

Proof: Let u be an eigenfunction of — A^w, with eigenvalue /i < v — ^. Then 
w = xnu G Dom(Ax~) satisfies (with ci = ^/j) 

This follows from exponential decay of w and is proved in the same way as (|82p . This 
implies |/i — A| < Ce^'^^'^ for some A S spec(— Ajjfoo) (in particular, ci may be re- 
placed by c). Since /i < v — ^, A must be an eigenvalue of — Ax°° • Now apply the 
Spectral Approximation Lemma [7] to A — — Aj^oo, W = spanjw}, with Ag = /i and 
e = Ce~'^^ , 6 = dist(/i, spec(— Ajfoo)\{A}). Since the interval (/i— (5, ^-l-(5) intersects 
the spectrum of — Ax~ only in A, we get dist(span{w}, Eig|;^} (— Ajfoo )) < Ce~'^^, 
and this implies (using exponential decay again) 

(84) dist(span{w}, Eig^PP ^) < Ce'^^. 
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Finally, applying this to an orthonormal basis of eigenf unctions of — A^^w with 
eigenvalues in /, we get (f83|) from (f37| . 

End of proof: The first statement of the theorem is contained in Step 2, and ((8T|) 
and dSSl) together imply ^ by Lemma [H:). □ 

7.2. Eigenvalues arising from the interior of the continuous spectrum. To 

define the approximate eigenspaces, recall Theorem If 51 Let 

= {af^AT : i = f , 2, . . . ; /o = f , . . . , Tj, a^jy G (0, amax]}, J^n = {i^+a'^ ■ a e ^w}- 

Recall that a E An iff HEa^^ satisfies the matching conditions for some f ^ 0. The 
corresponding space of ip is Va.N from (I65p . Therefore, the function 



(85) -B^^^ = 



on 

ni;^,^ + Xjvn^i^a.^ on 



is smooth on Xq and in the domain of Aj^n. 
For / C R let 

Eig^^w •= span{Ea,v, v + a^El, a E An, ^ E Pa.w}- 

Also, if / = (A, B) and 5 > then let 

h := (A-(5,B + ,5). 

We call dimT'o, TV the multiplicity of a resp. of A = + a^. 

Theorem 22. Assume — Ax°° ^as J^o -eigenvalues in [i^, A,„ax]- The numbers 
X E An are the approximate eigenvalues o/— A^w, with approximate eigenf unctions 

linear combinations of Ea,ip, V? S 'Pa,N ■ The errors are of order e^^^. 

More precisely, given a sufficiently small co > there are constants C, c > 
such that: 

a) Let fii be the eigenvalues bigger than iz+e^^"^"^ o/ — A^j-jv , arranged in increasing 
order and counted with multiplicity. Also, let Xi be the elements of An, arranged 
in increasing order and counted with multiplicity. Then, for all i, 

(86) |A,-M,| <Ce-=^. 

b) Let I CZ {ly + Ce~^^, Amax]- If there is no X E An in /2Ce-=" \ I then 

(87) dist3ymm(Eig,,^A„ Eig^PP) < Ce"^^, /' = Ice-^N . 

See Section 17.41 for the modifications needed in case there are L^-eigenvalues in 

[l^, Amax]- 

The statement in b) is complicated due to the possible crossings of the branches 
for different i in Figure [2] These do not occur on the line z = aN for bounded 
i (corresponding to fixed k as in Theorem [TJ, and one obtains: 

Corollary 23. The eigenvalues of —Axn form clusters of width Ce~'^^ around 
the X E An . For any Cq > there are c > 0, A^o such that for N > Nq the clusters 
around the X < v + CqN~'^ are disjoint and the span of eigenf unctions of — Aj^-n 
corresponding to the X-cluster has distance less than e^''^ from 

span{i?a^^ : ip E Via)} 

where X = v + , V{a) = ker(/ - e^°'^'^^aS{a)). 
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Proof. The first statement is just (j86|) . There are only a bounded number of A £ A^v 
satisfying A < + CqN^"^ by the Weyl asymptotics (l58|) . They are polynomially 
separated, i.e. there are X e N, c > (depending on Co) so that |A — A'| > cN~^ 
for any two different such values, by Theorem [15] and the fact that the are 
different analytic functions. This implies the disjointness of clusters for large N 
and that the separation condition in Theorem [22b) is satisfied for / = {A}, and 
this gives the last claim. □ 

Proof of Theorem \2^ As always, we write \ ^ v . 

Step 1: Show that the approximate eigenfunctions are actually such: For any 
intervals as in the theorem we have, for C sufficiently large, 

(88) dist(Eig^PP,Eig,,,v)<l- 

Proof: Let A € Km and let W\ — Eig^^j jy. For w = Ea,^ & W\ one has 
Aw = -Xw + 2L-^x'NdxIi^E^,^ + L-^'I^Il^Ea,,^, 

and since x'n^x'n ^^'^ supported in (N — 1,N), one obtains from ([?T|) . applied to 
u = Ea,^ with M € (A^ - 1, TV), that 

(89) ||(A^«+A)i.||;,« <Ce-^^||Hlx«- 

Since W\ C Dom(A_5fiv) we may apply the Spectral Approximation Lemma [7] to 
the operator A = — Aj^jv, with Aq = A and e = Ce~'^^ and J > to be chosen, and 
this gives 

(90) dist(Eig^PP ^, Eig(^„5^^+5),^) < Ce-=^/<5. 

If (5 = 206-"^ then this is ([88]) with / = {A}. 

To obtain (|88|) for arbitrary /, first observe that (I90p implies dist(Eig|Pj ^, Eigj jy) < 

Ce~^^/(5 for any interval J containing (A — 5, A + ^). Next, use Lemma [20] to- 
gether with a version of ([57]) for almost orthogonal subspaces to conclude that 
dist(EigfP,Eig^_, ^) < (dist(Eig^PP ,/5) + Ce-'=^). Since Ajv has at most 

0{N) elements, the left hand side is bounded by NCe"^^ / 5. Hence, choosing 
5 = 2A^Ce~^^ (or e"^ ^ with c' smaller than c and N large) one obtains ([88]l . 

Step 2: Show that each eigenvalue ijl> v + e^^^"^ of —/S.xg is exponentially close 
to some A g Ajv and that, under the assumptions of the theorem, 

(91) dist(Eig,,,jv, Eigfj^) < Ce--^. 

Proof: Let (A^^ « + = ^ = v + 0^ , so (3 > e'^""^ . For E = E^,^ recall the 
notation {E°,E^) = {\lE.:^=o,nd^E.:^=o) V y. V . ¥oy a An denote 

(92) Tc^^N = {{E\ E^): E = e Vo.,n)- 
Step 2a: {u^^u^) is close to Cp by the elliptic estimate: 

(93) dist(span{(w°, u^)}, Cp) < Ce""^. 

Proof: Apply the basic elliptic estimate. Lemma [T51 as follows: Let {E'^,E^) 
be the orthogonal projection of (m°,u^) to £p. This corresponds to a scattering 
solution E for spectral value /i. Let v — u ~ E. From (A + fi)E — 0, B±{fi)E = 
it follows that v satisfies ([74]) with f — 0, g = B±{fi)u, C'p — the orthogonal 
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complement of Cp and h = Q. From Lemma [19] it follows that ||5||Hi/2(y) < 
Ce^'^'^\\u\\H-i{x°)i and then (|75|) gives 

(94) |Mk2(xo) < Ce-^^||u||ff2(xa), v^u-E. 
This implies ||u||//2(j5fO) < C\\E\\fj2(^x°)i a-nd with (l76l) we get 

(95) ||u||h2(xo) <C||(i?o,Si)||vxy. 

Next, the trace theorem imphes w^)||yxv < C\\v\\u-2(^x°-^, and so (|M|) gives 

(96) |l(t;°,«i)|!yxv < Ce--''\\{E\E^)\\v.v. 

which after writing E = u — v and absorbing the v term into the left hand side 
gives \\{v^\v^)\\vy,v < Ce-=^||(u",Mi)||yxi', that is 1^. 

Step 2b: Use the matching conditions to show: ((93)l implies 

(97) dist(^,^jv) < Ce-'^^ 

(98) dist(span{(u",u^)}, J^c.tv) < Ce""^. 

aeAjv:|Q-/3|<Ce-=" 

The Ta,N are defined in ([9^ . 

Proof: Since {AxN+fi)v — one obtains from ([SG]) . using | cos /?Af| < 1, |/3^^ sin(/3iV)| < 
AT, that 

II MG[3iv)\\v < Uv'Wv + {N+ -)\\v')\\v < e=«^||(«",«i||yxv. 

a 

From MCf3{u) = we have MC/3(£') = - 1^10,3(1;), and then dH]) gives 

(99) II MGp{E)\\v = II MC^(t.)||y < Ce--'^\\{E",E^)\\v^v 

where ci = c — cq. If -E = Ep^^p then clearly ||(£'", i?^)||yxy < C||iy9||y, so we can 
apply the Stability Theorem [T7I with 6 = Ce~^i^ and obtain (I^T)) with c = ci, as 
well as dist(span{(^}, 0^, Vc^n) < Ce~^^^ where C2 = ci/(2 + 2dimF) and the 
sum is over a' e An satisfying |a — a'| < Ce~'^'-^^^. This implies that {E^, E^) and 
hence, by ([M]) . (it°,u^) has distance at most Ce^'^^^ from ©^/^q'.at, and hence 
(igg)) . with c = C2. 

End of proof of Step 2: The estimate ||w||x" 1^ C N\\{v'~' , v^)\\v xv for eigensolu- 
tions on which are a difference of a scattering solution and an eigenfunction on 
Xq (use (I5T|) and a modification of the derivation of (|75|) ) shows that ([M]) implies 

(100) dist({M}, Eig^PP) < Ce-'^, J = (a^ - Ce-^^, + Ce-=^). 

Finally, we apply this to an orthonormal set of eigenfunctions u with eigenvalues 
in r. Lemma II;) gives dist(Eig^,_^, Eig^PP^) < Ce'^^ with /" = (J')ce-" = 
^2Ce-<^" and hence (|91|) . since by assumption any A G /" n Ajv must already lie in 
/. 

End of proof of Theorem [22} and ([OTjl give part b) of the Theorem by 
Lemma [BJj). Part a) then follows easily: Since A^v has 0{N) elements, we may 
cover it by intervals Ik of length at most CNe~'^^ , satisfying the hypothesis of b) 
(note that any A G Ajy is at least ly + cN^"^ by (|55|) . (|64p since zi > 0). The /i^ 
must then be in the Ce~'^^-neighborhoods of the Ik by b), and this imphes a), with 
slightly smaller c. □ 
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7.3. Eigenvalues close to the threshold v. Recall from LemmalTlltliat ni?o, 
satisfies the matching conditions at a: = iff (T<i>+ = and = 0. Denote 

V]^ := V+ n F+, $1 = the projection to oi <^ e V. 

Note that, since a and iS'(O) do not commute in general, some care is needed with 
this notation. For example, usually ^ + (with <i>^ defined analogously). 
For $ e y+ let 

flOl) E .= /^o^*-" 

\ni;o,$,o + xwn^£^o,*.o on 

and 

(102) Eig^PP := {^o,*.o : * e K^}- 
Here we prove the following: 

Theorem 24. Suppose v is not an -eigenvalue of ~Ax°° ■ There are Cq > and 
c, C > such that all eigenvalues of —Ax^ in the interval {v — cq, + g-2co-/v^ ^^.g 
actually in I := {y — e~'^^ , v + e~'^^), and 

(103) dist„(Eig,_jv, Eig;^!'^) < Ce"=^. 

See Section 17. 4[ esp. (|123p . for the modification needed in case v is an L^- 
eigenvalue. 

Proof. Step 1: Show that the approximate eigenvalues are actually such: 

(104) dist„(Eig::pp,Eig,_^) < 1. 

Proof: This is proved in exactly the same way as ([88|) (with / = {A}). 

Step 2: Show that each eigenvalue /Lt e {v ~ co,u + e"^'^"^) of — A^n is in / and 
that 

(105) dist(Eig,^^, Eig^PP) < Ce-^^. 
Proof: Define 

(106) .fo,jv = {(i^°, E^): E = E„^^^„, $ e V+} 

analogous to (|92|) . Let u be an eigenfunction of — A^w, with eigenvalue fi. Let 

Since there are no scattering solutions with fj, < ly we compare u with a scattering 
solution for X ~ i>. 

Step 2a: (it°,M^) is close to Cq by the elliptic estimate: 

(107) dist(span{(u*',Mi)},£o) < S' := 5 + e''^ . 
Proof: Denote 

$ = u", vfr = u\ E ^ Eo,-S>+,^_- 
We apply the elliptic estimate, Lemma [THl to the difference v — u — E, with 
X = ly and £[, = {{v°,v^) : {v°)+ = 0,{v^)- = 0}. By this is transver- 

sal to £o- V satisfies (A + i')v = (A + u)u = {v — fi)u and Pco,£j| (w", f ^) = 
(since {v°,v^) € £'q by construction). Also, B±{iy)v = B±{i')u — B±{^)u + 
{Bi_{v) - B±{fi))u, and dTT]) gives ||B_L(Ai)-u||^i/2(y) < Ce"''^||u||H2(xo) while 
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clearly \\{Bj_{v) — _Bj^(/i))u||^i/2(y-) < C\iy ~ fj,\ ■ \\u\\h-^i^x°)- The elliptic estimate 
([75]) now gives 



-cN 



e 



(108) \\v\\hhx«)<CS'\\u\\h^^x'>), S':^S 
For S sufficiently small, this implies ||u||^2(x") < C\\E\\fj2(^xo}- Using (ffS)) and 
||(w°, 'i;^)||yxy < C||i^||ff2(xo) from the trace theorem we get IKw", i;^)||yxi^ < 
C5'\\{E°,E')\\vxv, that is, HMD. 

Step 2b: Use the matching conditions to show: p07|) implies that there is a 
constant co > so that for j.i > ly — cq 

(109) \fi - iy\ < Ce-"^ 

(110) dist(span{(u°,Mi)}, J^o.w) < Ce-"'^ . 

Proof: First, note that = T*_ (with T := |S"(0)), E^ = and 

1,0 = - £;0 = $_ - T^"., =u^ - E^ ^ imply that (fTUTl) is equivalent to 

(111) ||$_ - r*_|| + 1|*+|| < cs'i\\^+\\ + 

Using ($, ^f) = + one gets from this 

(112) !($,*)!< Cy||$||2 + C7||*||2. 

Note that this estimate does not involve the V± splitting. This is essential for the 
argument. 

We now consider the cases fi < v and fi > v separately. 
The case ^ < v 

For the sake of clarity we assume for the following argument that L = I. The 
case of general L requires only adjusting the constants. 

Let a = ^/S and t = tanh aiV. The matching conditions ([24|) for IIu = $ coshaa;+ 



(113) $"+a"4*"=0 

(114) ai<I>+ + 1'+ = 0. 

This implies ($,*) = ($+,*+) + = -at\\<!>+\\'^ - a-^ll^'-jp and so 

(115) K^-,*)! = at\\<^>+\\^ + a-H\\^-\\^ 

(116) ||$||2 = \\<^>+f + a-H^\\^-f 

(117) ||*|P=a2i2||^+||2_^||^-||2 

Now (|112l) implies that at least one of the following inequalities must hold: 

(118) at < CS' + Cah^ 

(119) a'H<CS'a-h^ + C. 

Multiply the second inequality by a^, plug in 6' = + e^^^ and use < i < 1 
to see that the second inequality implies the first. So (|118p holds. We claim that 
there is ao > so that a < qq implies a < Ce~'^^ . To see this, first observe that 
the Ca?t^ term on the right may be absorbed into the left for sufficiently small a, 
since t < 1. So we get at < Cc? + Ce~^^. Now for a > N^^ we have t > tanhl, 
so the term may be absorbed into the left, which yields a < Ce~'^'^ , while for 
a < we have t > c'aN for some constant c' > 0, and this gives a < Ce~^^^'^. 
We have shown that a^ ~ ijl\< Ce^"'^ if a < oq, that is, (|109p . 
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In particular, a ^ N. Now use (|112p again in conjunction with (I115p - (|117p . 
where we keep only the term on the left hand side, to obtain A^H^*"!]^ < 
Ce"'^^(||$+|P + A^2||^'"|P) + C||^'"||2. For large N the terms on the right may 
be absorbed, and one obtains ||^^|| < Ce~'^^||$+||. Together with (|113p this gives 

(120) ||$-|| < Ce-"^||$||. 
Also, since pli)) gives < Ce-'=^||$+||, we obtain 

(121) ||vf||<e-=^||$||. 
Now (imp implies 

(122) < e-"^||$||. 

Finally, (|120p and p22p imply by an elementary argument H^-^;};!! < Ce^'^^H^H. 
Therefore, ||($,^') - ($i,0)|| < Ce-''^\\{<i>,^)\\, that is, ITTOl) . 

The case > v. Let again a — v^, but now t = tanaTV, with a~^t := N 
if a = 0. Since we assume a < Ce~'^^ for this case, we may argue as in the 
last part of the argument for ji < v (starting with the paragraph before (|120p ). 
Observe that now (|114p is replaced by -ai$+ + = 0, which yields |($,*)| = 
1 0*11$+ IP — a~4||^'~||^| instead of ()115p (so pi2p gives only weaker conclusions 
than before), but the conclusions are still valid since atjl^+p < Ce~''^||<i>"'"|p. 

End of proof of Step 2: Exactly as in the proof of Theorem [^H it follows from 
Step 2b that dist(span{u}, Eigg^^^) < Ce~^^; applying this to an orthonormal basis 
of Eigj and using Lemma [5Jl),e) we get the claim. 

End of proof of Theorem I24t The claim follows directly from Steps 1 and 2, 
using Lemma [61d)- □ 

7.4. The case of embedded L^-eigenvalues. Here we sketch the modifications 
necessary in the arguments to deal with the case that — Ax°° has L^-eigenvalues 
embedded in the essential spectrum. For simplicity we will restrict to the analysis 
of eigenvalues near in case that v is an eigenvalue of — Ax°=. The case of Li^- 
eigenvalues bigger than v is treated similarly. 

Let n = {u& L^{X°°) : (Ax- + iy)u = 0}. If H ^ 0, Theorem [22 holds with 
the definition of EigJ^^^ replaced by 

(123) Eig;^PP {^o,$,o $ G V+} + {xnu -.ueH}. 

Also, Theorem [22] continues to hold as stated (if there are embedded eigenvalues 
X > ly then its statement has to be modified in a straightforward way). 

To prove this, we have to first modify the elliptic estimate. Lemma [THl We are 
interested in a near 0. Let Ho be the space of restrictions of elements of Ti to 
X° and Hq its orthogonal complement in L'^{X^). Then the elliptic estimate as 
stated cannot hold since the homogeneous problem (i.e., f = g = h = Oin (|74|) 1 
has solution space Ho- However, the same argument as given there shows that the 
same estimate holds if u G "Hq, and this gives 

(124) \\u - Pou\\h2(x'>) < Ci Wfh^xo) + ||.g||ffi/2(y) + \\h\\vxv ) 

where Pq : L^{X^) Ho denotes the orthogonal projection. 

Next, we have the following almost orthogonality statement analogous to Lemma 

m 
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If u is an eigenfunction of — Aj^-n , 



with eigenvalue n ^ v, then 



(125) 



I|-fou||ff2(xo) < C 




For the proof it suffices to show the same estimate for the L'^{X^) norm of Pqu by 



for all u' G TCq, with scalar product and norms in L^(X°). For this, do the analogous 
calculation as at the start of the proof of Lemma [20l then use that Ilu' = and 
that IIm^^jvIIi' — C'e^^^ljit'll , ||n-'-U3;=Ar||y < Ce^^^||w||, with the same estimate 
for the ^-derivatives. 

Now the proof of Theorem goes through as before since for /i > z/ + Ce~^"^ 
(jl25p shows that (|124p reduces to the 'old' elliptic estimate ([75)) . 

For the proof of Theorem [M] we first observe that Step 1 may be proved simply 
by a combination of the proofs of the Steps 1 in Theorems \n\ and Next, for 
Step 2 we may assume right away that |/i — < Ce~^^ since otherwise the 'old' 
elliptic estimate holds (see the previous paragraph) and the proof does not need to 
be modified. Now for an eigenfunction u of — Aj^-jv let mq € 7i be the eigenfunction 

of — Ajc°° restricting to Pqu, and let t2 = u — mq. Then, since ?i|xci = {I — Pq){u\x°) i 
(|124p is just the 'old' elliptic estimate for u, so the proofs of Steps 2a and 2b go 
through for u instead of u as before (with minor modifications when using (|77|) , and 
the matching conditions only satisfied up to an exponentially small error because 
of the mq term, which is inessential for the resulting estimate), and this gives that 
u is exponentially close to an -Eo,*,o and hence that u = uq + u \s exponentially 
close to Eigj^P^. 

7.5. Proof of Theorem[2j First, choose cq > so that the conclusion of Theorem 
[Ml holds. Since Vq = , (|103p gives the eigenvalues in Theorem [5]d), by Lemma 
O)) (and actually precise information on the eigenfunctions). Next, with this Cq 
apply Theorem [211 then ([SO)) gives the eigenvalues a) (for Tp < v), and Theorem 
[22l then ([86[) gives the eigenvalues in c). The eigenvalues close to those Tp which 
are > v are obtained using the argument in the preceding subsection. The cited 
theorems also give that there are no other eigenvalues. 



Proof of Theorem [^ We first discuss how to obtain the eigenvalues of a quantum 
graph. The metric graph (G, 2L) (that is, the graph G with given edge lengths 
2Ze, considered as a one-dimensional simplicial complex, i.e. as a union of intervals 
glued at the vertices) is just the space defined in with vertex and edge 

manifolds all equal to a point. Here we disregard the dimension requirement on 
the vertex and edge manifolds; but since the dimension requirement was never 
used (except implicitly in the validity of the theorems of scattering theory) we 
may use all previous results except those on scattering theory. Scattering theory 
is replaced as follows. A boundary condition at the vertices of G corresponds to a 
scattering matrix 5*0 (a), defined for a ^ by the requirement that the function 
e~'"^(pG -I- e^"'^ Sg{cx)^g on satisfy the boundary condition for each ipQ g Vq- 
By Lemma |12I this function satisfies the matching condition at a; = 1, i.e. extends 
to a smooth function on the metric graph, iff 



standard elliptic regularity, and for this we need to show | (m, u') \ < Cf- 



^NMI^'ll 



8. Identifying the quantum graph; special cases 



(126) 



(/-e 



ia2L 
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Since i^g = 0, this means that the positive eigenvalues of this quantum graph are the 
squares of those a for which ()126|) has a solution tpc (counted with multiplicity, 
defined as dimension of the space of those ifc)- 

On the other hand, from ([6]) and ((59)) we have that the positive bk in Theorem [T] 
are precisely the squares of those z > for which (/ — e*^^-^cr5(0))(^ has a solution 
</? 7^ 0, counted with multiplicity. 

It follows that we should choose boundary conditions for the quantum graph 
such that 

(127) Scia) = S{0) for aU a. 

In particular, we should take Vq — V, which leads us to consider functions which 
on the edge e take values in Afe', Lemma [TU] shows that (|127p yields the boundary 
conditions ([8]) , Finally, Lemma [Tl] shows that also the zero eigenvalues of the 
quantum graph correspond to the bk = 0- □ 

Note that the proof also gives a correspondence of the leading parts of eigen- 
functions (since they are determined by ipc and ip). 

We recover previously known results easily. The operator on the quantum graph 
in Theorem [3] is sometimes called the limit operator. For the following statement, 
see the remarks after that theorem, and for the notation the beginning of Section 

E 

Theorem 25. Suppose all vertex and edge manifolds are connected. Let Xq be the 
smallest eigenvalue of —Axo,^/, where M means that, in addition to the D/N (resp. 
Robin) boundary conditions at dX^ \ Y , we impose Neumann boundary conditions 
at Y . Then —Ax°° has no -eigenvalues < Xq, and: 

a) If X > V then we have Dirichlet conditions, i.e. decoupling, in the limit operator. 

b) If X = V = Q then we have Kirchhoff boundary conditions in the limit operator. 

In particular, for Neumann boundary conditions on all of OXq one has Kirchhoff 
boundary conditions, as proved in [5]. 

Proof. We prove the following stronger statement: If A < Aq then the equation 
(Ajfoo + X)u = can have no bounded solution, and if A = Aq = the only 
bounded solutions are constant on each X"^ . 

By Theorem [3l with Remark 3 following it, this implies the theorem since L^- 
solutions are bounded and since elements in the (+l)-eigenspace of 5(0) correspond 
to bounded solutions by Theorem [91;) . 

First, by Lemma El for a bounded solution u with A < wc must have if) = 
0, Vfe = OVA; in ^ resp. ^ and and this implies 

{ux=o, (9^w)x=o)y < 0. 
The same is true for an i^-solution for any A. Green's theorem implies 

(128) a / u(-AlZ) = -(M,=o,(a«u).=o)y + / iVwp 

Jxo Jxo Jxo 

(129) > / |Vu|2, 

Jxo 

SO < A if u\xo ^ 0. Since may be taken as test function in the 

variational characterization of Ao, this implies Ao < A and hence the first claim. If 
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A = Aq = then it implies Vu = 0, so u is constant on (since X° is connected) 
and hence on X!^ by unique continuation. Note that for = connectedness of Ye 
imphes that A^. = C canonicaUy (the constant functions). □ 

9. Appendix: Monotone unitary families 

In this appendix we coUect some resuhs on analytic one-parameter families of 
unitary operators which we need. Discussion and proofs can be found in 

Let U{x) be a family of unitary operators on a Hermitian vector space V, of 
dimension AI < oo, depending real analytically on a: G R. Then 

(130) D(x) -U'(x)U(xy^ 

I 

is symmetric, where U'{x) is the derivative with respect to x. Assume that U is 
monotone, i.e. D{x) is positive for all x, and more precisely that there are constants 

dmin, rfmax, ^2 > SUCh that 

(131) dnuni < D{x) < dmax/, \\U" {x)\\ < for aU X. 
Denote 

W{x) = Ker(/ - U{x)) and 2 = {x : W^x) ^ {0}}. 

Thus X e Z iff U{x) has eigenvalue one. 

A special case of this setup is U{x) = e"J7o for a unitary Uq. Then Z is discrete 
and 27r-periodic, and W{x) is the eigenspace of Uq with eigenvalue e~*^. The 
following statements generalize this and well-known facts about eigenspaces to our 
more general situation. 

Lemma 26. Z dM. is a discrete subset, and more precisely for all A < B 

rB 



(132) 



^ dimW^(a;)-— / trD{x)dx 



< M(:= dhnV) 



x:A<x<B 

The following lemma mimics the independence of the eigenspaces. 

Lemma 27. Let I he an interval of length at most '^f'Jj ■ Then the spaces W{x), 
X Cz I, are independent, i.e. 

(133) // (fix €E W{x) for each x ^ I r\ Z and Lp^ — ^ then = ^ Vx. 

X 

The following lemma gives a stable version of almost orthogonality. 
Lemma 28. Assume ip \ Q satisfies 

(134) ||(/-C/(a;o))^||<£||^||. 
Then 

(135) dist(xo,Z)< '^^ 



dmin 

Furthermore, there is a constant C only depending on dmim dmax, d2, M such that 
if s < C^^ then, with Pw denoting the projection to ®\x-xo\<^/£^^^^' 

(136) ||^-Pv^^||<CeW/+i)||^||. 
We also need a fact about 2-parameter perturbations. 
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Theorem 29. Let U (x, y) he a unitary operator in a finite- dimensional Hermitian 
vector space depending real analytically on x,y G M. Assume 

(137) 7^^"^^° a<(xo,yo). 

Then the set {{x,y) : U{x,y) has eigenvalue one} is, in a neighborhood of {xo,yo), 
a union of real analytic curves x = Xj{y). The corresponding projections Pj{y) 
to the eigenspace of U{xj{y),y) with eigenvalue one are also analytic functions 
of y ^ yo, extending analytically to y = yo, and '^jPjiyo) is the projection to 
keT{I -U{xo,yo))- 
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